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We propose that the topological magneto-electric (ME) effect, a hallmark of topological insulators 
(TIs), can be realized in thin films of TIs in the v = 0 quantum Hall state under magnetic field 
or by doping two magnetic ions with opposite signs of exchange coupling. These setups have the 
advantage compared to previously proposed setups that a uniform configuration of magnetic held or 
magnetization is sufficient for the realization of the topological ME effect. To verify our proposal, 
we numerically calculate ME response of TI thin films in the cylinder geometry and that of effective 
2D models of surface Dirac fermions. The ME response is shown to converge to the quantized value 
corresponding to the axion angle 6 = ±7r in the limit of the large top and bottom surface area of 
TI films, where non-topological contributions from the bulk and the side surface are negligible. 


I. INTRODUCTION 

Topological aspects of the condensed matter systems 
have been a subject of intensive studies since the dis¬ 
covery of the quantum Hall effect (QHE) [1, 2]. The 
topological nature of the QHE is encoded in the Hall 
conductance, which is exactly quantized in units of e 2 //i 
to the Chern number characterizing the global topology 
of wave functions constituting a manifold in the Hilbert 
space. The Hall current supporting the quantized Hall 
conductance flows through the edge channels along the 
boundary of the sample (bulk-edge correspondence). 

The extension of the QHE to time-reversal symmetric 
systems is realized in topological insulators (TIs) [3, 4]. 
Here, the Z 2 topological index characterizes the global 
topology of wave functions of Kramers pairs. Helical edge 
channels and surface Dirac fermions exist on the bound¬ 
ary of two-dimensional (2D) and three-dimensional (3D) 
TIs, respectively, as a consequence of nontrivial topology. 
The spin Hall conductance of 2D TIs is an analog of the 
Hall conductance in the QHE. However, the spin current 
is not conserved in the presence of spin-orbit interactions 
which are indispensable for TIs; the spin Hall conduc¬ 
tance is neither a well-defined nor quantized quantity. 
The charge conductance along the helical edge channels 
in 2D TIs was shown to be quantized at 2 e 2 /h [5], but the 
quantization is not so accurate as in the QHE. Further¬ 
more, one cannot relate the charge conductance to the 
Z 2 topological index in 3D TIs. The quantized quantity 
characteristic of TIs, corresponding to the Hall conduc¬ 
tance in the QHE, has been looked for. 

The topological magneto-electric (TME) effect [6-15] 
has been suggested as a phenomenon that is directly 
related to the Z 2 topological index for 3D TIs. The 
TME effect is a quantized response of polarization to 
applied magnetic fields or a quantized response of mag¬ 
netization to applied electric fields [6, 7]. The quantized 
polarization current was proposed to be accessible by 
transport measurements [6-8]. Alternatively, the TME 
effect was proposed to be observed by optical measure¬ 


ments of Faraday or Kerr rotation showing a quantized 
response [10, 11]. However, no experimental observation 
of the TME effect has been reported so far. 

The TME effect of a 3D TI is theoretically described 
by the axion action (with c = 1), 

S = S~J = ^ / ■ B, 

(1) 

where the angle 0 is equivalent to the Z 2 index: 9 = 
7r (mod 27 t) for TIs and 0 = 0 (mod 27 t) for trivial insu¬ 
lators [6]. Transforming Eq. (1) to the integral over the 
surface of a 3D TI yields the Chern-Simons action, 

Op 2 1' 

Scs = 8^nj dtd2x ^ XA ^ A N ( 2 ) 

which describes the QHE on the TI surface with the Hall 
conductance a xy = (0/2ir)(e 2 /h). For TIs with 9 = ir, 
this amounts to a xy = e 2 /2 h. Microscopically, the half¬ 
integer QHE ( <j xy = ±e 2 /2 h) on the TI surface occurs 
once surface Dirac fermions acquire a mass and have the 
Hamiltonian 

H = vf{-VxO v + PyCT x ) + M<Jz, (3) 

where vf is the Fermi velocity. The mass gap M of the 
surface Dirac fermions can be introduced, for example, 
through the exchange coupling J of Dirac fermions to 
the magnetization m (pointing perpendicular to the sur¬ 
face) of doped ferromagnetic ions; M = Jm. In order 
to realize the TME effect described by Eqs. (1) and (2), 
two conditions must be satisfied, (i) The magnetization 
points outward or inward over the whole surface of a 3D 
TI such that surface Dirac fermions are fully gapped on 
any surface, (ii) The Fermi energy must be tuned to be 
inside the gap induced by the magnetization to. This 
magnetic configuration is very difficult to realize experi¬ 
mentally since the external magnetic field favors the same 
magnetization direction for both top and bottom sur¬ 
faces. In addition, the fine tuning of the Fermi energy is 
still very difficult experimentally although some groups 
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FIG. 1. Setups for the TME effect in (a) TI thin films in the 
v = 0 QHE, and (b) magnetically doped TI thin films. Right 
panels are schematics of energy bands of surface states. 


have succeeded [16, 17]. We note in this regard that lo¬ 
calization effect can relax the second condition. That 
is, when the surface Dirac fermions are localized by mag¬ 
netic impurities, the Fermi energy can be set in the wider 
energy window of localized states [14]. Furthermore, the 
magneto-electric cooling effect can help the magnetiza¬ 
tion m align perpendicular to any surface because of the 
energy gain proportional to the volume which eventually 
overcome the surface energy [14]. However, there is no 
experimental study to pursue this possibility up to now. 

In this paper we propose alternative routes to realizing 
the TME effect in thin films of 3D TIs by either applying 
external magnetic field or doping two kinds of magnetic 
impurities. The first approach utilizes the QHE of TI 
thin films [18-20]. In particular, we are interested in the 
novel v = 0 QHE which was recently observed in the pres¬ 
ence of a finite potential difference between the top and 
bottom surfaces [20]. In this case, the edge channels are 
gapped on the side surfaces due to the finite thickness 
effect, which can produce interesting spintronics func¬ 
tions [21], We point out that this v = 0 quantum Hall 
system provides an ideal laboratory to realize the TME 
effect, as discussed in detail below. We also propose that 
the TME effect can be achieved by doping Cr and Mn to 
upper and lower halves of TI thin films, respectively. 


II. TOPOLOGICAL MAGNETO-ELECTRIC 
EFFECTS IN TI THIN FILMS: BASIC IDEA 


The TME effect described by the actions in Eqs. (1) 
and (2) predicts that charge polarization P is induced by 


applying a magnetic field B, 


P = 


A 

2-7T h ’ 


( 4 ) 


where 9 = tt (mod 2-7t). In this section we explain how 
this quantized response of charge polarization is obtained 
for the proposed two setups of TI thin films by discussing 
charge response of quantum Hall states of Dirac fermions 
at the top and bottom surfaces. 

Let us consider a thin film of a TI in the magnetic field 
B = (0,0, B) which is applied perpendicular to the film; 
see Fig. 1(a). In this case, Dirac fermions at the top and 
bottom surfaces form Landau levels (LLs). For simplicity 
we assume that the Fermi energy is at the Dirac point, 
Ef = 0. In this case, important low-energy excitations 
are states in the LLs at the Fermi energy on the top and 
bottom surfaces, while the excitations are gapped in the 
bulk and at the side surface (the latter due to the finite¬ 
thickness effect). In addition, we introduce a potential 
energy difference between the top surface (Uq > 0) and 
bottom surface (—Uq). For small magnitude of the po¬ 
tential asymmetry Uq , the LLs at the top and bottom 
surfaces are occupied up to Nt = — 1 and Nb = 0 LLs, 
respectively, as indicated in Fig. 1(a), where Nt and Nb 
are LL indices on the top and bottom surfaces. This sit¬ 
uation realizes the v = 0 QHE which was experimentally 
observed in Ref. [20]. More generally, the v = 0 QHE is 
realized in TI thin films when LLs at the top and bot¬ 
tom surfaces are filled up to Nt = — N — 1 and Nb = N 
LLs with any integer N, where N can be controlled by 
changing the potential asymmetry Uq. Then the charge 
densities at the top and bottom surfaces (tit and ns) 
become 


" T = ( w+ 5)x |iS| - " ii = -( Ar+ 5)7r |B| - (5) 

If the top and bottom surfaces are located at z = d and 
z = —d, respectively, the charge polarization along the z 
direction is given by 

p = ^[d nT + (-d) nB }= (n+^J^\ b \- (6) 

The polarization in the v = 0 QHE coincides with the 
value expected from the TME effect, Eq. (4) with 8 = 
(2 N + l)7r, for B > 0. For B < 0, the axion angle is 
given by 9 = —(2N + l)-7r. This jump of 9 at B = 0 
takes place because the surface bands become gapless at 
B = 0. In both cases, the axion angle takes a nontrivial 
value 9 = tt modulo 27 t. We emphasize again that, in 
the setup with the v = 0 QHE, an arbitrary odd integer 
9/ 7T = 2 N + 1 can be realized by tuning the potential 
asymmetry Uq with gating the top and bottom surfaces 
independently such that the highest occupied LLs are set 
to Nt = — N — 1 and Nb = N. 

Next we consider the second setup in which the upper 
and lower halves of a TI film are doped with two types 
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of magnetic ions, e.g., Cr and Mn. Here a crucial as¬ 
sumption is that the two kinds of magnetic ions have ex¬ 
change couplings J of opposite signs with Dirac fermions 
in the TI film. This is indeed the case with Cr and Mn, 
as can be seen from the fact that Cr-doped and Mn- 
doped Bi 2 Te 3 -based TIs show quantized anomalous Hall 
conductivity of opposite signs when the Fermi energy is 
tuned within the excitation gap of surface Dirac fermions 
[16, 17, 22]. Then a uniform magnetization depicted in 
Fig. 1(b) produces an exchange field pointing outward or 
inward at the top and bottom surfaces, and the desired 
constant Dirac mass Jma z in Eq. (3) is achieved except 
for the side surface. In the 3D geometry in Fig. 1(b), 
massive Dirac fermions show quantum anomalous Hall 
effect a xy = — sgn( Jm)e 2 /2h and a xy = sgn( Jm)e 2 /2h 
at the top and bottom surfaces, respectively. Since the 
electron density on each surface is given by n = a xy B , 
applying a magnetic field induces the charge polarization 


p =r/ d 


—sgn(Jra) — 


B 




= —sgn (Jm) — B. 


(7) 


Thus the TME effect with the axion angle 9 = 

— 7T sgn( Jm) is realized with a uniform magnetization. 
Since Dirac fermions on the side surface are gapped be¬ 
cause of the finite thickness, the TME effect should be ob¬ 
served when the chemical potential is tuned into the en¬ 
ergy gap of the side-surface excitations. Compared with 
the original proposal with the non-uniform configuration 
of magnetization [14] which is experimentally difficult to 
realize, this setup has the advantage that it only requires 
a uniform magnetization, which facilitates experimental 
realization. 

For both setups the charge polarization can be detected 
as a charge current in response to time-dependent B. 


III. NON-TOPOLOGICAL CONTRIBUTIONS 

FROM THE BULK AND THE SIDE SURFACE 

In general, the axion angle can deviate from 9 = 0,7r 
(mod 27 t) when both time reversal symmetry T and inver¬ 
sion symmetry P are broken. The v = 0 QHE is achieved 
in TIs by applying magnetic field, which breaks T, and 
by inducing potential asymmetry between the top and 
bottom surfaces, which breaks P. Thus the axion angle 
in the v = 0 QHE setup is not necessarily quantized. 
In the previous section we have shown that polarization 
corresponding to the quantized value of 6 = 7r (mod 27r) 
is obtained when we neglect contributions from the bulk 
and focus on contributions from the top and bottom sur¬ 
faces. While the charge polarization in the TME effect 
is accompanied by a flow of the surface current, non- 
topological contributions may also arise from the current 
flowing through the bulk. Furthermore, the side surfaces 
connecting the top and bottom surfaces in Fig. 1 may 


give rise to non-topological contributions. Then the nat¬ 
ural question is whether charge polarization of TI thin 
films is really quantized as Eq. (4) predicts with the ax¬ 
ion angle 9 = ir (mod 27r) when contributions from the 
bulk and the side surface are taken into account. To ad¬ 
dress this question, we consider the 3D bulk of TI thin 
films and calculate the charge polarization and the cur¬ 
rent distribution induced by increasing B. 

Suppose that the TI sample has the cylinder geometry 
shown in Fig. 2, in which the top and bottom surfaces 
are located at z = ±d and the side surface is at r = R 0 . 
In order to describe the bulk electronic states of the TI 
film, we take a simple massive Dirac Hamiltonian 

H = H 0 + V, Hq = Vf(p ■ cr )t x + m(r )r~, (8) 

where p = ( p x ,p y ,p z ) is the momentum, and cr = 
(cr x , a y , a z ) is the spin Pauli matrices, while t x . t z are 
the orbital Pauli matrices. We assume that the chemical 
potential p = 0. The Dirac mass (i.e., the bulk band gap) 
m(r) = rriin > 0 inside the TI film. The Dirac mass m(r) 
changes its sign at the surface, and m(r) = m ou t < 0 out¬ 
side the TI. In the setup of the v = 0 QHE in Fig. 1(a), 
V is given by 


V(r) = vj?eA ■ (tt x + U(z), (9) 

where A is the vector potential in the symmetric gauge 
A = ^B(—y,x, 0), and the potential energy U(z) = Uq 
for 2 > 0 and U(z) = —Uq for z < 0. In the setup of a 
TI doped with Cr and Mn in Fig. 1(b), V is given by 

V(r) = M(z)a z , (10) 

where M(z) = M 0 for z > 0 and M(z) = —M 0 for z < 0. 
In the cylindrical coordinates (r, tp, z), wave functions are 
written as 

ip(r, p, z) = r - 2 e l (^ _ -r)v , ^( r . ) z ) (11) 

for the angular momentum t £ Z + i. Then the Dirac 
equation is reduced to [H 0 (£) + V](/)e = E(£)<j>e, where 

H 0 (£) =v F ^p r a x + y <7 y +p z a z ' S jr x +m(r)T z , (12) 

and <j)(r, z ) obeys the Dirichlet boundary condition <j> = 0 
at r = 0. 

Small magnetic field 5B(t) applied in the z direction 
(in addition to the constant magnetic field B for the v = 0 
QHE case) induces a change in the polarization along the 
z direction that is proportional to SB(t). The change 
in the polarization leads to polarization current J pro¬ 
portional to ddB(t)/dt. From the linear response theory 
(see Appendix for details), the induced current density is 
written as 

J( r o) ST 2h Im K^b'( T, o)l^j)(V , jl c> l^)] 

dSB/dt ^ {Ei - Ej) 2 ’ 1 1 

j eu 
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FIG. 2. Current density (indicated by arrows) induced by 
small time-dependent magnetic field J(r) / (dSB/dt) for (a) 
TIs in the v = 0 QHE with B = 10T, Uo = 20meV, and 
(b) magnetically doped TIs with Mo = lOmeV. Insets are 
schematic pictures of the cylinder geometry of TI films. We 
adopted parameters: vf = 5 x 10 5 m/s, 2d = 10 nm, Ro = 
40 nm, min = 250 meV, and m ou t = —750 meV. 
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FIG. 3. Schematic picture of the effective 2D model with 
coordinates (X, Y) along the surface of a TI thin film. 


current in the z direction flows to connect current distri¬ 
butions at the top and bottom surfaces, and it is exactly 
the polarization current along the z direction in the TME 
effect. From the discussions focusing on the top and bot¬ 
tom surfaces in the previous section, the polarization cur¬ 
rent is expected to be equal to the time derivative of P 
in Eq. (4) with the axion angle at 9 = n for the v = 0 
QHE from Eq. (6) with N = 0, and at 9 = —tt for the 
magnetically doped TI from Eq. (7) with Jm = Mq > 0. 
However, the axion angles obtained from the integration 
of the polarization current plotted in Fig. 2, 

S = ? -.RUMB/it) J = °>> < 14 > 

are found to be 0.4 and —0.3 for the v = 0 QHE and the 
magnetically doped TIs, respectively, and both deviate 
from the quantized value 9/2t: = ± 7 ). We attribute these 
deviations to finite-size effects of the surface of TI films, 
considering that the bulk contribution to the polarization 
current is negligibly small as shown in Fig. 2(a) and (b). 
Specifically, Dirac fermions on the side surface give non- 
topological contributions to the polarization current be¬ 
cause of broken time-reversal and inversion symmetries. 


where i and j label occupied states ipi (i € O) and 
unoccupied states (j £ U) with energies Ei and 

Ej, respectively, j(r 0 ) = — ev F S(r — r 0 )crT x , and O = 
\ev F (-ycr x + x<j v )t x . 

The current density in Fig. 2 shows that the cur¬ 
rent flows along the surface of the TI film, not through 
the bulk. Thus the charge polarization in the present 
setup is indeed a topological phenomenon carried by 
the Dirac fermions on the surface of TIs. The surface 
current of topological nature is understood as follows. 
Time-dependent SB induces an electric field along the 
azimuthal (ip) direction. Then the Hall current flows 
along the r-direction on the top and bottom surfaces. 
The directions of the Hall current are opposite for the 
top and bottom surfaces, where a xy = ±e 2 / 2 h and 
a xy = Te 2 /2 h, respectively. On the side surface, the 


IV. ANALYSIS FROM SURFACE WEYL 
FERMIONS AND THE QUANTIZED AXION 
ANGLE 

The current distribution in Fig. 2 indicates that the 
polarization current is negligible in the gapped bulk and 
flows only at the surface of TI thin films. We now focus on 
low-energy Dirac fermions at the surface and study their 
contribution to the magneto-electric effect. In particular, 
we are interested in finite-size corrections which cause 
the deviation of 9 from the quantized value. For this 
purpose, we introduce an effective 2D model (Fig. 3) for 
surface Dirac fermions, in which surfaces of a 3D TI thin 
film are extended on the 2D (A, Y) plane. The whole 
2D surface has linear dimensions of (2 L x + Sd)L y and 
the y direction is assumed to be periodic. The original 
3D coordinates (x, y , z) of each surface are related to 2D 
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coordinates (X, Y) as follows: 

(—X + d + Tpj Y, —d), 
(X-d-^,Y,d), 
{x,y,z) = \ (-%SY,X), 


(a) 


(^,Y,L x + d-X), 


f,Y,-L x -d-X), 


—L x - d < X < —d, 
d -X .£/.£ + d, 

|X| <d, 

+ d < X, 

X < - d, 

(15) 


see Fig. 3. 

For the v = 0 QHE setup the effective Hamiltonian for 
the surface Dirac fermions is given by 

H = vf { —p x & y + [py + eA Y (X)]a x } + U(X)a 0 . (16) 

Here, we take the Landau gauge Ay{X) = Bx(X). The 
potential difference between the top and bottom surfaces 
is given by the potential 

Uq , d < X < L x + d, 

U(X) = { -U 0 , —L x — d < X < -d, (17) 

0, otherwise, 

as schematically illustrated in Fig. 4(a). Since the 2D 
model is translationally invariant along the Y direction 
under the Landau gauge, the wave number ky{= k y ) is 
a good quantum number. 

The polarization in the z direction is given by 


p =2srjw- 


(18) 


We note that the spectrum of z is bounded in our 2D ef¬ 
fective model so that the above equation is well defined. 
From the linear-response theory, the axion angle 9 deter¬ 
mined from the charge polarization 9/2n = hSP/e 2 SB is 
given by 
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2d 



9 —e h f dk y 2 

27r 2 dL x e 2 J 2n ^ Ei, ky - E ] k 

j eu 

x Re[('0j > fc H \z\ipj,k y )(ipj,ky \evFXO’ x \il)i'k y )\, (19) 

where i and j label occupied states i/)i,ky (* € O) and 
unoccupied states ipj,k y (j € U) with momentum k y and 
energy E i ky and Ej t k y , respectively. To compute Eq. 
(19) numerically, we discretize the 2D model. In order 
to avoid the fermion doubling on the lattice, we used a 
regularization scheme in the momentum space. Namely, 
we used the plane wave basis for the X direction labeled 
by kx and set a cutoff to the wave number kx as \kx\ < 
A. The numerical results are obtained by extrapolating 
the cutoff A —>• oo. 

We show results for the axion angle computed from Eq. 
(19) in Fig. 4. The axion angle 9/2tt approaches 1/2 as 
L x — y oo. This shows that non-topological contributions 
from the side surface become negligible as the ratio L x /d 
is increased. Indeed the deviation of the computed axion 


FIG. 4. (a) Schematic picture of a TI thin film under magnetic 
fields, (b) Band structure of a TI thin film for L x = 200nm. 
(c) Axion angle in the v = 0 QHE plotted against the system 
size L x . The inset shows axion angle plotted against the in¬ 
verse of the system size 1 /L x . We performed the calculation 
with parameters: vf = 5 x 10 5 m/s, B = 10T, Uo = 20meV, 
and 2d = lOnm. 


angle 9/2n from the quantized value 1/2 scales with the 
ratio of the areas of the side surface and the top/bottom 
surfaces, d/L x , as shown in Fig. 4(c). Thus the axion 
angle converges to the quantized value 1/2 in the limit 
where the top/bottom surfaces have much larger area 
than side surfaces. Since typical values of the thickness 
2d and the size L x of TI films are 2d = 10 nm and L x = 
1 mm [20], the axion angle 9/2n is expected to show a 
quantized behavior in practice. 

Next we move on to the 2D surface model for the mag- 
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netically doped TI. The effective Hamiltonian for the sur¬ 
face Weyl fermions is given by 

H = V F {-p X <7y + PYCTx ) + M ■ <T, (20) 

where we assumed no potential asymmetry between the 
top and bottom surfaces. The upper and lower halves 
of the TI film are doped with Cr and Mn, respectively. 
A uniform magnetization along the z direction induces 
effective exchange fields +M 0 and — M 0 along the 2 di¬ 
rection in the upper and lower halves of the TI film, re¬ 
spectively, because signs of exchange couplings are op¬ 
posite for Cr and Mn. In the 2D model of surface Weyl 
fermions, these exchange fields are perpendicular to the 
top and bottom surfaces 

M(X) = (0,0, Mq) (21) 

for | A" ± (L x / 2 + d)| < L x j 2 and produce a mass term 
Mqc z , while they are perpendicular to the side surface 

M(X) = (±M 0 ,0,0) (22) 

for —d < ±X + L x + d < 0 and 0 < ±X < d and enter as 
a vector potential ±Mocr x in Eq. (20). This is schemati¬ 
cally illustrated in Fig. 5(a). We numerically obtain the 
polarization and the axion angle using the linear response 
theory [Eq. (19)] in the same manner as in the v = 0 QHE 
case. The result is shown in Fig. 5. We again see asymp¬ 
totic behavior 0/2t: —> —1/2 with L x —> 00 for fixed d. 
In a similar way to the case of the v = 0 QHE, the devi¬ 
ation of the axion angle from the quantized value —1/2 
scales with d/L x as shown in Fig. 5(d). Thus there is a 
non-topological contribution proportional to the area of 
the side surface, but it diminishes in the thermodynamic 
limit L x /d —> 00 where typical experiments of thin films 
are performed. 


V. DISCUSSIONS 




ky (nm' 1 ) 
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We have proposed that the TME effect should be ob¬ 
served in TI films in the v = 0 QHE and in magnet¬ 
ically doped TI films. We have shown that the polar¬ 
ization current in these setups are topological in that 
bulk contributions are negligible and show a quantized 
response in the limit of large top/bottom surface area. 
The non-topological contribution to the polarization cur¬ 
rent is proportional to the area of the side surface, which 
is interpreted as the axion angle 6 at the side surface de¬ 
viating from the quantized value ±7r due to the breaking 
of time-reversal and inversion symmetries. For a thin- 
film sample of area 1mm 2 and a H-field sweeping speed 
of IT/s, the polarization current is estimated to be 20pA. 
This magnitude of the polarization current is accessible in 
the current state of transport measurements [23]. Appli¬ 
cation of ac magnetic fields induces a polarization current 
oscillating in time in the 2 -direction. This ac polarization 
current can be detected, for example, with electrodes ca- 
pacitively coupled to top and bottom surfaces and mea¬ 
suring currents between two electrodes as schematically 


FIG. 5. (a) Schematic picture of a TI thin film doped with 
magnetic ions, Cr and Mn. (b) Band structure of a TI thin 
film for L x = 200nm. (c) Axion angle plotted against the 
system size L x . The inset shows axion angle plotted against 
the inverse of the system size 1/L X . We performed the calcu¬ 
lation with parameters: vf = 5 x 10 5 m/s, Mo = lOmeV, and 
2d = lOnm. 


shown in Fig. 6. Besides the corrections to the TME ef¬ 
fect [Eq. (4)] due to broken T and P symmetries in our 
setup, there can be other corrections from other mech¬ 
anisms, such as the coupling between top and bottom 
surfaces. However, the dominant part of the polarization 
should come from the topological contribution. 

In the ideal situation where the axion action in Eq. (1) 
with 0 = 7r is realized over the whole sample, the total 
charge induced by B is proportional to the cross sec¬ 
tional area of the sample perpendicular to B as follows. 
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FIG. 6. A setup to observe the TME effect in TI thin films 
by detecting the polarization current induced by changing 
magnetic fields. 


We suppose that B is parallel to the 2 direction, and 
the sample geometry is defined by height functions for 
upper and lower surfaces, z = z u (x,y) and z = zi(x,y ), 
respectively, for x, y within the projected region S of the 
sample onto the xy-plane. The total charge on the upper 
surface is given by 


1 


j-z u (x,y) q g 2 

Q= I dxdy— -- - -- / dz——B 

J s Zu{x,y) - zi(x,y) J Zl ( x , y ) 2tt h 


9 e 2 

~ 2W i B S ’ 


(23) 


where As = f s dxdy is the cross sectional area of the 
sample in the xy-plane. Conversely, the axion angle 9 
deviates from the quantized value in our setups for the 
TME effect, since there is non-topological contribution 
proportional to the area of the side surface as shown in 
Fig. 4(c) and Fig. 5(c). 

Note added: After the completion of this work, a paper 
on a related topic by Wang et al. [24] has appeared. 
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Appendix: Linear response theory for current 
density and polarization 

We give a brief derivation of Eqs. (13) and (19) from 
linear response theory. The Hamiltonian 


value of an operator B in the ground state of H , 

(B) = B„e~™, (A.2) 

which is expanded in powers of A. The Fourier compo¬ 
nent B u in linear order in A is obtained from the linear 
response theory, 


B u = —, 


dr e luJT tr 


e - iHoT / h [A, p 0 }e iHoT/h B 


(A.3) 


where the density matrix po = |0){0| is the projection to 
the many-body ground state |0) of H 0 . For noninteract¬ 
ing fermions, po can be written as po = ILeo l'0*)(^*l 
in terms of the occupied single-particle states ipi. The 
Fourier component B u in Eq. (A.3) is then rewritten as 


B u = — 


dr 


X [ 

i e o 
jeu 


K . A ..pi{ E i-Ej+tiu)T/h 

A-* 2 7 *'■'11'-' 


—Bn A, 


’jlSXljl 


i (Ej — E i + huj ) r / h 


= E 


ieO 

jeu 


BijAji 


_ | (BjjAji)* 

Ei — En + Hoj Ei — Ej — Tux 


- E 


2 [(Ei - Ej)Re(BijAji) - ihwlm(BijAji )] 


i e o 
jeu 


(Ei - Ejf - (huf 


(A.4) 

where O and U stand for sets of occupied and unoccupied 
single-particle states in the ground state of Hq, and Aij 
and Bij are matrix elements of the operators A and B. 
We have used the hermiticity of the operators A and B 
when passing from the first line to the second line. In 
the dc limit Eq. (A.4) reduces to 


B : 


'w-VO 


= E 


i e o 
jeu 


2R e(BijAji) 
Ei — Ei 


(A.5) 


The axion angle 6 is defined in Sec. IV as a response 
coefficient of polarization P to magnetic field B. The 
linear response formula [Eq. (19)] is obtained by setting 
B = P and A = evFXa x B in Eq. (A.5). 

In order to obtain the formula in Eq. (13) for the cur¬ 
rent density that is proportional to the time derivative of 
the magnetic field, we need to extract from (B) the com¬ 
ponent that is proportional to the time derivative of the 
perturbation, dA(t)/dt = — iuiAe~ lult . That component 
is given by 


lim R e[B u /(—iuj)\ 

CD—>-0 


2hlm( Bij Aj i ) 

,r 0 ^-^> 2 

jeu 


(A.6) 


H = H 0 + A(t), (A.l) 

consists of an unperturbed part Hq and a perturbation 
A(t) = Ae~ lult . We are interested in the expectation 


When the magnetic field is varied as B —> B + SB , the 3D 
bulk Hamiltonian H in Eq. (8) is perturbed by 5B O = 
5 evF(—yo’ x + x<j v )t x 5B. Equation (13) is obtained by 
setting B = j and A = O in Eq. (A.6). 
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